Simultaneous self-supervised reconstruction and denoising of sub-sampled MRI data with Noisier2Noise

Most existing methods for Magnetic Resonance Imaging (MRI) reconstruction with deep learning assume that a high signal-to-noise ratio (SNR), fully sampled sampled dataset exists and use fully supervised training. In many circumstances, however, such a dataset does not exist and may be highly impractical to acquire. Recently, a number of self-supervised methods for MR reconstruction have been proposed, which require a training dataset with sub-sampled k-space data only. However, existing methods do not denoise sampled data, so are only applicable in the high SNR regime. In this work, we propose a method based on Noisier2Noise and Self-Supervised Learning via Data Undersampling (SSDU) that trains a network to reconstruct clean images from sub-sampled, noisy training data. To our knowledge, our approach is the first that simultaneously denoises and reconstructs images in an entirely self-supervised manner. Our method is applicable to any network architecture, has a strong mathematical basis, and is straight-forward to implement. We evaluate our method on the multi-coil fastMRI brain dataset and find that it performs competitively with a network trained on clean, fully sampled data and substantially improves over methods that do not explicitly remove measurement noise.Comment: Submitted to IEEE International Symposium on Biomedical Imaging (ISBI) 202

PEAR: PEriodic And fixed Rank separation for fast fMRI

In functional MRI (fMRI), faster acquisition via undersampling of data can improve the spatial-temporal resolution trade-off and increase statistical robustness through increased degrees-of-freedom. High quality reconstruction of fMRI data from undersampled measurements requires proper modeling of the data. We present an fMRI reconstruction approach based on modeling the fMRI signal as a sum of periodic and fixed rank components, for improved reconstruction from undersampled measurements. We decompose the fMRI signal into a component which a has fixed rank and a component consisting of a sum of periodic signals which is sparse in the temporal Fourier domain. Data reconstruction is performed by solving a constrained problem that enforces a fixed, moderate rank on one of the components, and a limited number of temporal frequencies on the other. Our approach is coined PEAR - PEriodic And fixed Rank separation for fast fMRI. Experimental results include purely synthetic simulation, a simulation with real timecourses and retrospective undersampling of a real fMRI dataset. Evaluation was performed both quantitatively and visually versus ground truth, comparing PEAR to two additional recent methods for fMRI reconstruction from undersampled measurements. Results demonstrate PEAR's improvement in estimating the timecourses and activation maps versus the methods compared against at acceleration ratios of R=8,16 (for simulated data) and R=6.66,10 (for real data). PEAR results in reconstruction with higher fidelity than when using a fixed-rank based model or a conventional Low-rank+Sparse algorithm. We have shown that splitting the functional information between the components leads to better modeling of fMRI, over state-of-the-art methods

Adapting water management to climate change in the Murray–Darling Basin, Australia

Climate change is threatening water security in water-scarce regions across the world, challenging water management policy in terms of how best to adapt. Transformative new approaches have been proposed, but management policies remain largely the same in many instances, and there are claims that good current management practice is well adapted. This paper takes the case of the Murray–Darling Basin, Australia, where management policies are highly sophisticated and have been through a recent transformation in order to critically review how well adapted the basin’s management is to climate change. This paper synthesizes published data, recent literature, and water plans in order to evaluate the outcomes of water management policy. It identifies several limitations and inequities that could emerge in the context of climate change and, through synthesis of the broader climate adaptation literature, proposes solutions that can be implemented when basin management is formally reviewed in 2026

Probabilistic optimization for conceptual rainfall-runoff models: a comparison of the shuffled complex evolution and simulated annealing algorithms

Automatic optimization algorithms are used routinely to calibrate conceptual rainfall-runoff (CRR) models. The goal of calibration is to estimate a feasible and unique (global) set of parameter estimates that best fit the observed runoff data. Most if not all optimization algorithms have difficulty in locating the global optimum because of response surfaces that contain multiple local optima with regions of attraction of differing size, discontinuities, and long ridges and valleys. Extensive research has been undertaken to develop efficient and robust global optimization algorithms over the last 10 years. This study compares the performance of two probabilistic global optimization methods: the shuffled complex evolution algorithm SCE-UA, and the three-phase simulated annealing algorithm SA-SX. Both algorithms are used to calibrate two parameter sets of a modified version of Boughtoh's [1984] SFB model using data from two Australian catchments that have low and high runoff yields. For the reduced, well-identified parameter set the algorithms have a similar efficiency for the low-yielding catchment, but SCE-UA is almost twice as robust. Although the robustness of the algorithms is similar for the high-yielding catchment, SCE-UA is six times more efficient than SA-SX. When fitting the full parameter set the performance of SA-SX deteriorated markedly for both catchments. These results indicated that SCE-UA's use of multiple complexes and shuffling provided a more effective search of the parameter space than SA-SX's single simplex with stochastic step acceptance criterion, especially when the level of parameterization is increased. Examination of the response surface for the low-yielding catchment revealed some reasons why SCE-UA outperformed SA-SX and why probabilistic optimization algorithms can experience difficulty in locating the global optimum.Mark Thyer and George Kuczera, Bryson C. Bate

Time-encoded pseudo-continuous arterial spin labeling: Increasing SNR in ASL dynamic angiography

Purpose: Dynamic angiography using arterial spin labeling (ASL) can provide detailed hemodynamic information. However, the long time-resolved readouts require small flip angles to preserve ASL signal for later timepoints, limiting SNR. By using time-encoded ASL to generate temporal information, the readout can be shortened. Here, the SNR improvements from using larger flip angles, made possible by the shorter readout, are quantitatively investigated. Methods: The SNR of a conventional protocol with nine Look-Locker readouts and a 4 (Formula presented.) 3 time-encoded protocol with three Look-Locker readouts (giving nine matched timepoints) were compared using simulations and in vivo data. Both protocols were compared using readouts with constant flip angles (CFAs) and variable flip angles (VFAs), where the VFA scheme was designed to produce a consistent ASL signal across readouts. Optimization of the background suppression to minimize physiological noise across readouts was also explored. Results: The time-encoded protocol increased in vivo SNR by 103% and 96% when using CFAs or VFAs, respectively. Use of VFAs improved SNR compared with CFAs by 25% and 21% for the conventional and time-encoded protocols, respectively. The VFA scheme also removed signal discontinuities in the time-encoded data. Preliminary data suggest that optimizing the background suppression could improve in vivo SNR by a further 16%. Conclusions: Time encoding can be used to generate additional temporal information in ASL angiography. This enables the use of larger flip angles, which can double the SNR compared with a non-time-encoded protocol. The shortened time-encoded readout can also lead to improved background suppression, reducing physiological noise and further improving SNR

Time-encoded pseudo-continuous arterial spin labeling: Increasing SNR in ASL dynamic angiography

Purpose: Dynamic angiography using arterial spin labeling (ASL) can provide detailed hemodynamic information. However, the long time-resolved readouts require small flip angles to preserve ASL signal for later timepoints, limiting SNR. By using time-encoded ASL to generate temporal information, the readout can be shortened. Here, the SNR improvements from using larger flip angles, made possible by the shorter readout, are quantitatively investigated. Methods: The SNR of a conventional protocol with nine Look-Locker readouts and a 4 (Formula presented.) 3 time-encoded protocol with three Look-Locker readouts (giving nine matched timepoints) were compared using simulations and in vivo data. Both protocols were compared using readouts with constant flip angles (CFAs) and variable flip angles (VFAs), where the VFA scheme was designed to produce a consistent ASL signal across readouts. Optimization of the background suppression to minimize physiological noise across readouts was also explored. Results: The time-encoded protocol increased in vivo SNR by 103% and 96% when using CFAs or VFAs, respectively. Use of VFAs improved SNR compared with CFAs by 25% and 21% for the conventional and time-encoded protocols, respectively. The VFA scheme also removed signal discontinuities in the time-encoded data. Preliminary data suggest that optimizing the background suppression could improve in vivo SNR by a further 16%. Conclusions: Time encoding can be used to generate additional temporal information in ASL angiography. This enables the use of larger flip angles, which can double the SNR compared with a non-time-encoded protocol. The shortened time-encoded readout can also lead to improved background suppression, reducing physiological noise and further improving SNR

Development and Validation of the Gene Expression Predictor of High-grade Serous Ovarian Carcinoma Molecular SubTYPE (PrOTYPE).

PURPOSE: Gene expression-based molecular subtypes of high-grade serous tubo-ovarian cancer (HGSOC), demonstrated across multiple studies, may provide improved stratification for molecularly targeted trials. However, evaluation of clinical utility has been hindered by nonstandardized methods, which are not applicable in a clinical setting. We sought to generate a clinical grade minimal gene set assay for classification of individual tumor specimens into HGSOC subtypes and confirm previously reported subtype-associated features. EXPERIMENTAL DESIGN: Adopting two independent approaches, we derived and internally validated algorithms for subtype prediction using published gene expression data from 1,650 tumors. We applied resulting models to NanoString data on 3,829 HGSOCs from the Ovarian Tumor Tissue Analysis consortium. We further developed, confirmed, and validated a reduced, minimal gene set predictor, with methods suitable for a single-patient setting. RESULTS: Gene expression data were used to derive the predictor of high-grade serous ovarian carcinoma molecular subtype (PrOTYPE) assay. We established a de facto standard as a consensus of two parallel approaches. PrOTYPE subtypes are significantly associated with age, stage, residual disease, tumor-infiltrating lymphocytes, and outcome. The locked-down clinical grade PrOTYPE test includes a model with 55 genes that predicted gene expression subtype with >95% accuracy that was maintained in all analytic and biological validations. CONCLUSIONS: We validated the PrOTYPE assay following the Institute of Medicine guidelines for the development of omics-based tests. This fully defined and locked-down clinical grade assay will enable trial design with molecular subtype stratification and allow for objective assessment of the predictive value of HGSOC molecular subtypes in precision medicine applications.See related commentary by McMullen et al., p. 5271.Core funding for this project was provided by the National Institutes of Health (R01-CA172404, PI: S.J. Ramus; and R01-CA168758, PIs: J.A. Doherty and M.A.Rossing), the Canadian Institutes for Health Research (Proof-of-Principle I program, PIs: D.G.Huntsman and M.S. Anglesio), the United States Department of Defense Ovarian Cancer Research Program (OC110433, PI: D.D. Bowtell). A. Talhouk is funded through a Michael Smith Foundation for Health Research Scholar Award. M.S. Anglesio is funded through a Michael Smith Foundation for Health Research Scholar Award and the Janet D. Cottrelle Foundation Scholars program managed by the BC Cancer Foundation. J. George was partially supported by the NIH/National Cancer Institute award number P30CA034196. C. Wang was a Career Enhancement Awardee of the Mayo Clinic SPORE in Ovarian Cancer (P50 CA136393). D.G. Huntsman receives support from the Dr. Chew Wei Memorial Professorship in Gynecologic Oncology, and the Canada Research Chairs program (Research Chair in Molecular and Genomic Pathology). M. Widschwendter receives funding from the European Union’s Horizon 2020 European Research Council Programme, H2020 BRCA-ERC under Grant Agreement No. 742432 as well as the charity, The Eve Appeal (https://eveappeal.org.uk/), and support of the National Institute for Health Research (NIHR) and the University College London Hospitals (UCLH) Biomedical Research Centre. G.E. Konecny is supported by the Miriam and Sheldon Adelson Medical Research Foundation. B.Y. Karlan is funded by the American Cancer Society Early Detection Professorship (SIOP-06-258-01-COUN) and the National Center for Advancing Translational Sciences (NCATS), Grant UL1TR000124. H.R. Harris is 20 supported by the NIH/National Cancer Institute award number K22 CA193860. OVCARE (including the VAN study) receives support through the BC Cancer Foundation and The VGH+UBC Hospital Foundation (authors AT, BG, DGH, and MSA). The AOV study is supported by the Canadian Institutes of Health Research (MOP86727). The Gynaecological Oncology Biobank at Westmead, a member of the Australasian Biospecimen Network-Oncology group, was funded by the National Health and Medical Research Council Enabling Grants ID 310670 & ID 628903 and the Cancer Institute NSW Grants ID 12/RIG/1-17 & 15/RIG/1-16. The Australian Ovarian Cancer Study Group was supported by the U.S. Army Medical Research and Materiel Command under DAMD17-01-1-0729, The Cancer Council Victoria, Queensland Cancer Fund, The Cancer Council New South Wales, The Cancer Council South Australia, The Cancer Council Tasmania and The Cancer Foundation of Western Australia (Multi-State Applications 191, 211 and 182) and the National Health and Medical Research Council of Australia (NHMRC; ID199600; ID400413 and ID400281). BriTROC-1 was funded by Ovarian Cancer Action (to IAM and JDB, grant number 006) and supported by Cancer Research UK (grant numbers A15973, A15601, A18072, A17197, A19274 and A19694) and the National Institute for Health Research Cambridge and Imperial Biomedical Research Centres. Samples from the Mayo Clinic were collected and provided with support of P50 CA136393 (E.L.G., G.L.K, S.H.K, M.E.S.)

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer.

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC

Polygenic Risk Modelling for Prediction of Epithelial Ovarian Cancer Risk

Funder: Funding details are provided in the Supplementary MaterialAbstractPolygenic risk scores (PRS) for epithelial ovarian cancer (EOC) have the potential to improve risk stratification. Joint estimation of Single Nucleotide Polymorphism (SNP) effects in models could improve predictive performance over standard approaches of PRS construction. Here, we implemented computationally-efficient, penalized, logistic regression models (lasso, elastic net, stepwise) to individual level genotype data and a Bayesian framework with continuous shrinkage, “select and shrink for summary statistics” (S4), to summary level data for epithelial non-mucinous ovarian cancer risk prediction. We developed the models in a dataset consisting of 23,564 non-mucinous EOC cases and 40,138 controls participating in the Ovarian Cancer Association Consortium (OCAC) and validated the best models in three populations of different ancestries: prospective data from 198,101 women of European ancestry; 7,669 women of East Asian ancestry; 1,072 women of African ancestry, and in 18,915 BRCA1 and 12,337 BRCA2 pathogenic variant carriers of European ancestry. In the external validation data, the model with the strongest association for non-mucinous EOC risk derived from the OCAC model development data was the S4 model (27,240 SNPs) with odds ratios (OR) of 1.38(95%CI:1.28–1.48,AUC:0.588) per unit standard deviation, in women of European ancestry; 1.14(95%CI:1.08–1.19,AUC:0.538) in women of East Asian ancestry; 1.38(95%CI:1.21-1.58,AUC:0.593) in women of African ancestry; hazard ratios of 1.37(95%CI:1.30–1.44,AUC:0.592) in BRCA1 pathogenic variant carriers and 1.51(95%CI:1.36-1.67,AUC:0.624) in BRCA2 pathogenic variant carriers. Incorporation of the S4 PRS in risk prediction models for ovarian cancer may have clinical utility in ovarian cancer prevention programs.</jats:p